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Abstract 

I present a model in which dark energy interacts with matter. The former is represented 
by a variable equation of state. It is shown that the phantom crossing takes place at zero 
redshift, moreover, stable scaling solution of the Friedmann equations is obtained. I show 
that dark energy is most probably be either generalized phantom energy or the generalized 
Chaplygin gas, while phantom energy is ruled out as a dark energy candidate. 
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1 Introduction 



Numerous cosmological observations show that our universe is pervaded with a mysterious 
dark energy which constitutes more then seventy percent of the cosmic energy density p] . 
The presence of dark energy is essential to get a spatially flat universe as is observed today. 
Although dark energy is dominant in the universe, however we don't know its composition 
and origin. Theoretically there are several candidates that can account for dark energy 
including the cosmological constant, quintessence, Chaplygin gas and phantom energy etc, 
to name a few. However no single candidate among these can completely and satisfactorily 
meet the observational data. Another problem associated with it is that the exact form 
of the equation of state (EoS) of dark energy is yet to be discovered. In literature, 
several exotic EoS are proposed (see for example [2]). We in this paper will discuss the 
dynamics of dark energy with a recently proposed EoS that combines three forms of dark 
energy equations of state including generalized Chaplygin gas, phantom energy and its 
generalized form. Thus this single model presents unified cosmological model of analyzing 
dark energy dynamics. 

From the cosmological observations, it is now clear that there are three major cosmic 
ingredients including dark matter, radiation and dark energy. It is naturally expected that 
these cosmic species may entertain some form of interaction at some stage in the evolution 
of the universe. This interaction is specifically important in the present cosmic setting 
when we observe that the energy densities of dark matter and dark energy are almost of 
the same orders of magnitude. Thus it excludes radiation from the interacting system. 
Recently there have been several cosmological models build up having dark energy of the 
forms mentioned above interacting with matter [3]. It is interesting to note that most 
of these models predict matter dominated scenario in the late evolution of the universe 
where dark energy decays into matter. This conclusion apparently avoids the imminent 
cosmic singularities like the big rip or the big crunch. Thus it naturally predicts a cyclic 
form of the evolution of the universe from matter domination to dark energy domination, 
back and forth. 

In the present paper, I model my system by first assuming the spatially flat background 
spacetime. Further, I assume the spacetime to contain only dark matter and dark energy 
and having exotic interaction. Both dark matter and dark energy are specified by their 
respective equations of state. I choose a specific form of the interaction term that is 
well motivated from the dimensional considerations. The exact form of the interaction 
term will be originally motivated from the theory of quantum gravity which is as yet 
in the process of construction by a number of theorists. That term contains a coupling 
parameter whose numerical value determines the strength of the interaction while its 
signature describes the direction of flow of the energy. I shall assume its numerical value 
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to be between and 1 while its signature to be positive (showing flow of energy from 
dark energy to matter). Finally I perform stability analysis of the dynamical system and 
discuss my results. 



2 Modeling of dynamical system 

I start by assuming the background to be spatially flat, homogeneous and isotropic 
Friedmann- Robert son- Walker (FRW) spacetime 

ds 2 = -dt 2 + a 2 (t)[dr 2 + r 2 (d6 2 + sin 2 9 d(j> 2 )}. (1) 

The equations of motion corresponding to FRW spacetime filled with the two component 
fluid (matter and dark energy) are 

K 2 

H = - — (Pde + Pde + Pm), (2) 
O 

H 2 = ^(Pde + Pm). (3) 

Here k 2 = 8nG is the Einstein's gravitational constant and H = H(t) = a/a is the Hubble 
parameter. Also subscripts 'm' and 'de' correspond to matter and dark energy respectively 
while p and p are respectively pressure and energy density. In Eq. (2), matter is assumed 
to be pressureless (p m = 0) while the dark energy's equation of state is considered of the 
form [1] 

p de = -Br 2 ^p a de , a^O, (4) 

where B is a positive constant and a can be either positive or negative. Note that a = 1 
represents phantom energy, a > gives generalized phantom energy while a < repre- 
sents generalized Chaplygin gas. Thus the generalized phantom energy and Chaplygin 
gas models can be unified by the above equation of state. 

The energy conservation equation for the dynamical system under consideration is 

Pde + Pm + 3H (p m + p d e + Pde) = 0. (5) 

Due to interaction between the two components, the energy conservation would not hold 
for the individual components, therefore the above conservation equation will break into 
two equations: 

Pde + 3H(p de + Pde) = ~Q, (6) 

p m + 3Hp m = Q. (7) 
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Here Q is the interaction term which in general is a function of energy densities and the 
Hubble parameter i.e. Q(Hp m , Hpd e ) which upon Taylor expansion yields the form given 
below in (8). This interaction is demanded to be positive to avoid the violation of the 
laws of thermodynamics [5J. I also insert a dimensionless coupling parameter c in Q to 
determine the strength of the interaction, thus I have 

Q = 3Hc(p de + 

To study the dynamics of the system, I proceed by setting 

x = ]na = -ln(l + z), (9) 

which is termed as the e-folding time parameter and z is the redshift parameter. Moreover, 
the density and pressure of dark energy can be expressed by dimensionless variables u 
and v as 

n Pde ^Pde K 2 p de 

"^^^-p-^VP' V= MP- (10) 

The EoS parameter ujde is conventionally defined by 

Pde / 1 1 \ 

UJde = , (11) 

Pde 

which after using (10) becomes 

Ude = -■ (12) 

u 

The density parameters of dark energy and dark matter are related as 

n m = ^ = i-n de = i- u . (13) 

Using Eqs. (8)-(13) in (6) and (7), I obtain the following autonomous system 

du 

— = -3c-3v + 3uv, 14 
ax 



dv / v c 

— = _3 av (i + _ + _ 

dx V u u 



)+3«(« + l). (15) 



Note that the above dynamical system matches with [7] if a < 0. Their it is shown that 
interacting GCG can explain late time acceleration. I am here particularly interested 
in the cases when a = 1 and a > 0. The stability of this system is performed by first 
finding its critical points i.e. the points for which the left hand sides of the system vanish 
identically. Thus I obtain the only critical point: 

u c = 1 - c, (16) 
v c = -1. (17) 
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It needs to be checked whether the system is stable about (u c , v c ) . For this purpose, I 
linearize the system by substituting u = u c + 6u, and v — v c + Sv, where 5u and Sv are 
small perturbations about the critical point. I finally obtain 



dSu 



3v c Su + 3(-l + u c )Sv, (18) 
3 

— [a(Vc + cv c )5u + [(1 + 2v c — a)u 2 c — 2au c v c — acu c ]5v}. (19) 



dx 
d5v 



dx 



c 



Now the given dynamical system is stable about the critical point if the real parts of 
the eigenvalues are both negative. In this case, the critical point is called a stable node. 
However, if one is positive then such a point is termed a saddle point. If both positive, it 
is called an unstable node. The eigenvalues are 



It can be shown easily that the real parts of both the eigenvalues are negative. In Figures 
1 and 2, I show phase space diagrams of the solutions of my dynamical system for different 
choices of initial conditions. Here I took a = 0.07 and a = 1, representing generalized 
phantom energy and the phantom energy, respectively. Note that I am not discussing the 
case of a < as it already has been discussed in [7j. In figures 3 and 4, I provide pictorial 
representation of the solutions against the e-folding time parameter. It is interesting 
to note that the solutions are stable (without fluctuations or oscillations of any kind). 
Figure 5 shows an unstable solution which is consequently of not much physical interest. 
However it puts an upper limit on the choice of c = 0.6 if the dark energy is of phantom 
type a = 1. In Figures 6 and 7, I plot the dark energy state parameter against e-folding 
time parameter x. The former one presents a smooth transition from the quintessence 
to phantom regime while this transition takes place at x — 0, which corresponds to the 
present epoch. The later figure shows a rather stale case where the state parameter 
remains unchanged throughout the variation of x. I also point out here that this last 
result is unaffected for various values of c. Hence it can be said that dark energy can 
most probably be either generalized Chaplygin gas or the generalized phantom energy. 
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(20) 
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3 Conclusion 



Models of interacting dark energy possess tremendous potential to alleviate several cosmo- 
logical puzzles, most notably the cosmic coincidence problem and the phantom crossing 
scenario. The model predicts a dynamical nature of dark energy where its equation of 
state and energy density both vary over cosmic time which turns out to be consistent with 
the observations. As seen in the previous section, the phantom crossing scenario is well 
explained if the coupling parameter is non-zero but very small positive number compared 
to unity. In this paper, I have analyzed the model of interacting dark energy using an 
interesting equation of state which was originally specified in the context of wormhole 
physics, very similar to the Chaplygin gas that was originally motivated from an engi- 
neering problem. The equation of state amalgamates three different forms of dark energy 
and therefore its inclusion in the interacting model makes the model more comprehensive 
and unified. From the analysis of this paper, it can be concluded that dark energy could 
most probably of the form of generalized phantom energy or generalized Chaplygin gas 
while phantom energy is ruled out. Also, this work serves as a generalization of earlier 
studies in this direction [6], U\ E] which dealt with Chaplygin gas only. 
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Figure 1: The phase diagram of the interacting dark energy model with c = 0.7 and a = 0.07. 
The curves correspond to the initial conditions u(— 2) = 1.2, v(— 2) = —0.2 (green); u(— 2) = 
1.3,v(-2) = -0.3 (blue); u(-2) = 1.4,«(-2) = -0.4 (red); n(-2) = 1.5,«(-2) = -0.5 (brown). 
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Figure 2: The phase diagram of the interacting dark energy model model with c 
a = 1. The curves correspond to the initial conditions as given in Fig. 1. 
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Figure 3: The parametric functions u and v are plotted against the e- folding time parameter x. 
The parameters are fixed as c = 0.7 and a = 0.07. The initial condition is u(0) = 2, v{0) = —2. 
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Figure 4: The parametric functions u and v are plotted against the e-folding time parameter x. 
The parameters are fixed as c = 0.06, a = 1. The initial condition is u(0) = 2,v(0) = —2. 
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Figure 5: The parametric functions u and v are plotted against the e- folding time parameter x. 
The parameters are fixed as c = 0.6, a = 1. The initial condition is u(0) = 2,v(0) = —2. 
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Figure 6: The EoS parameter ujx = v/u is plotted against x. The parameters are fixed as c = 0.7 
and a = 0.07. The initial conditions are u(0) = 1, v(0) = —1 (solid line); n(0) = 1.1, v(0) = —1.1 
(dots); n(0) = 1.2, w(0) = -1.2 (dashes); n(0) = 1.3, u(0) = -1.3 (dash dot). 



13 



-0.8- 



v/u -1.0- 



-1.2- 



-1.4- 

~U 1 , 1 , 1 , 1 , 1— 

-0.4 -0.2 0.0 0.2 0.4 

x 



Figure 7: The EoS parameter ujx = v/u is plotted against x. The parameters are fixed as 
c = 0.5 and a = 1. The initial conditions are the same as in Fig. 6 
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